This is my PM566 Final Project website.


Introduction

Vaccine hesitancy, which is the reluctance or refusal to vaccinate despite the availability of vaccines, has been proposed by the World Health Organization (WHO) as one of the ten threats to global health in 2019 https://www.who.int/news-room/spotlight/ten-threats-to-global-health-in-2019, even before the global pandemic of COVID-19.

Vaccine is currently one of the most cost-effective ways of preventing disease outbreak and reducing disease consequence, and consequently disease’s impact on economic. As of 2019, 2-3 million deaths a year are prevented by vaccination. If global coverage of vaccinations is improved , an estimated 1.5 million of death could be further avoided.

As for the reason of vaccines hesitancy are complex, could be a result of numerous factors including lack of education, knowledge of how vacines are produce, distrust of government or authorities, or even just fear of needles. Its existence dates back the time vaccines were invented, however, it is only recently that it started to resurface to the discussion of mainstream media.

After the global outbreak of COVID-19 and the invention of COVID-19 vaccines, vaccine hesitancy plays an even more important role than before. We can observe vaccines being not just a measure of disease prevention, but also a mean of politics, yet we can also see a subset of population refused to get vaccinated even with the resources at hand.

This report is not to discuss the cause of vaccine hesitancy or the legitimacy of the claims by certain “anti-vaccinationism”, but to see if vaccine hesitancy actually impacts the death and case number of COVID-19 on a state level in the US.

Methods

Datasets

A total of 3 datasets will be used in this project. The first two are from US Center for Disease Control and Prevention, and the third one are from United States Census Bureau (an official website of the United State Government).

  1. The fisrt dataset is “Vaccine Hesitancy for COVID-19: County and local estimates” from https://data.cdc.gov/Vaccinations/Vaccine-Hesitancy-for-COVID-19-County-and-local-es/q9mh-h2tw

  2. The second dataset is “United States COVID-19 Cases and Deaths by State over Time” from https://data.cdc.gov/Case-Surveillance/United-States-COVID-19-Cases-and-Deaths-by-State-o/9mfq-cb36

  3. And for the comparison between states summary statistics to be reasonable and effective, I need the population data of different places. “County Population Totals: 2020-2021” from https://www.census.gov/data/tables/time-series/demo/popest/2020s-counties-total.html

Estimated state vaccine hesitancy level

The hesitancy level

Results

Tab 1

Tab 2

Run data

source("process_COVID_data.R")
p1_scatter <- cv_states_today %>% 
  plot_ly(x = ~pop_density, y = ~deathsper100k,
          type = 'scatter', mode = 'markers', color = ~state,
          size = ~population, sizes = c(5, 70), marker = list(sizemode='diameter', opacity=0.5),
          hoverinfo = 'text',
          text = ~paste( paste(state, ":", sep=""), paste(" Cases per 100k: ", per100k, sep="") , paste(" Deaths per 100k: ",
                        deathsper100k, sep=""), sep = "<br>")) %>%
  layout(title = "Population-normalized COVID-19 deaths vs. population density",
                  yaxis = list(title = "Deaths per 100k"), xaxis = list(title = "Population Density"),
         hovermode = "compare")

# filter out "District of Columbia"
cv_states_today_scatter <- cv_states_today %>% filter(state!="District of Columbia")

p2_scatter <- cv_states_today_scatter %>% 
  plot_ly(x = ~pop_density, y = ~deathsper100k,
          type = 'scatter', mode = 'markers', color = ~state,
          size = ~population, sizes = c(5, 70), marker = list(sizemode='diameter', opacity=0.5),
          hoverinfo = 'text',
          text = ~paste( paste(state, ":", sep=""), paste(" Cases per 100k: ", per100k, sep="") , paste(" Deaths per 100k: ",
                        deathsper100k, sep=""), sep = "<br>")) %>%
  layout(title = "Population-normalized COVID-19 deaths vs. population density",
                  yaxis = list(title = "Deaths per 100k"), xaxis = list(title = "Population Density"),
         hovermode = "compare")

Showcasing

Tab 1

Tab 2

Copyright © 2022, Tzu Yu Huang.